INVERSE EIGENVALUE PROBLEMS FOR PERTURBED 
SPHERICAL SCHRODINGER OPERATORS 
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Abstract. We investigate the eigenvalues of perturbed spherical Schrodinger 
operators under the assumption that the perturbation q(x) satisfies xq(x) £ 
L 1 (0, 1). We show that the square roots of eigenvalues are given by the square 
roots of the unperturbed eigenvalues up to an decaying error depending on the 
behavior of q(x) near x = 0. Furthermore, we provide sets of spectral data 
which uniquely determine q(x). 



1. Introduction 

Given a Schrodinger operator with a rotationally symmetric potential, separation 
of variables leads to the spherical Schrodinger operator (e.g., [TS], [H]) 

(1-1) H = -f^ + l -^±^ + q(x), 1 = 0,1,2,.... 

ax ar 

In this note we are interested in the case where the particle is confined to a finite ball 
(of radius 1 for notational simplicity — which can always be achieved by scaling x). 
This problem has attracted much interest in the past and several results concerning 
the eigenvalues of these problems have been derived. The first results are based 
on the seminal work by Guillot and Ralston [BJ which deals with the case I = 1 
and q £ L 2 (0, 1). It was later extended by Carlson [4], [5] and recently by Scrier 
[14] who extended their results to arbitrary 1 e No. However, the assumption 
q £ L 2 (0,1) clearly excludes the physically interesting case of a Coulomb type 
singularity q{x) = + ■ ■ ■ ■ This case was included for I = in the work of 
Savchuk and Shkalikov [T^], [13] who considered q £ W~ 1:2 (0, 1). Their work was 
later extended by Albeverio, Hryniv, and Mykytyuk who first covered the case 
q £ W~ 1 ' p (0, l),p6 [1, oo), for I = in [2] and later on extended this to all I e N 
in [3] using the double commutation method. Their condition includes in particular 
the case xq(x) € L p (0, 1), p £ [l,oo), which will be the condition we are interested 
in here. 

Our main motivation for this paper is the paper by Zhornitskaya and Serov [18 
who treat the general case I > — \ under the assumption q € i 1 (0, 1). We want 
to extend their results in several ways: First of all we will replace the condition 
q(x) e L 1 ^, 1) by xq{x) £ L l (0, 1). Moreover, they show in [18] that the Dirichlet 
eigenvalues satisfy 

(1-2) Hn = (ii+i ,„ + £n) 
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where i; + i )n = 7r(n + |) + 0(n 1 ) are the zeros of the Bessel function J i+ i/ 2 (z) 
and the error satisfies \s n \ < j (this is claimed for all n but only proven for large 
n). We will show that the error satisfies 

(1.3) e n = o(f^ d y). 

(For I — — i one has to replace q(y) by (1 — log(y))q(y).) In particular, for q £ 
i 1 (0, 1) we get e n = 0(n _1 ) and for a Coulomb type singularity q(x) — ^+i 1 (0, 1) 
we get e n = 0(rt _1 log(n)). 

Based on this information we will give some sets of spectral data which, in 
addition to the Dirichlet spectrum, uniquely determine q and I again generalizing 
the corresponding result from [18] to the case xq{x) £ L 1 (0,1). Moreover, in 
their construction they use the fact that a certain Wronskian does not vanish. 
Unfortunately this Wronskian can indeed vanish (we will give a simple counter 
example due to Ralston in Remark I2.7[) but we will show that the use of this fact 
can be avoided. 



2. The spherical Schrodinger operator 

Our prototypical example will be the spherical Schrodinger operator given by 

d 2 1(1 + 1) , s 1 

(2.1) Hl = -— + ±—±, xe(0,l), l>--. 



Note that we explicitly allow non-integer values of I such that we also cover the 
case of arbitrary space dimension n > 2, where 1(1 + 1) has to be replaced by 
1(1 + n - 2) + (n- l)(n - 3)/4 QS Sec. 17.F]. 

With the usual boundary conditions at x = (for I £ [—5,5)) and x = 1 

(2.2) \im x l ((l + l)f(x)-xf'(x))=0, f(l) = or /'(l) + 0/(1) = 0, 

x— >0 

it gives rise to a self-adjoint operator in the Hilbert space L 2 (0, 1). Two linearly 
independent solutions of 

-f"{x)+ l -^-f(x) = zf(x) 



(pi(z,x) = z 2i 4 +1 J^-J l+ i(^ x ), 



W*x) ~ \ log(z) J l+ 1 (v^a;), i + 

where Ji+1/2 and Y/+1/2 are the usual Bessel and Neumann functions. All branch 
cuts are chosen along the negative real axis unless explicitly stated otherwise. If / 
is an integer they of course reduce to spherical Bessel and Neumann functions and 
can be expressed in terms of trigonometric functions (cf. e.g. [TJ [17] and also [15j 
Sect. 10.4]). 
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Using the power series for the Bessel and Neumann functions one verifies that 
they have the form 

(2-6) <l> l {z,x)=x l + 1 T{l ^l )2l+l f l {zx% 
(2.7) 



Oi(z,x) 



r(z + f)2' +1 



2iTT9i(« 2 ), Z + ±eK+\N , 

[(log(2) - 7 k(^ 2 ) - log(x)/ ; (zx 2 ), I = -|, 

where fi(z), gi(z) are entire functions with /;(0) = gi(0) = 1 and 7 is the Euler- 
Mascheroni constant. 

In particular, both functions are entire and according to [TJ (9.1.16)] their Wron- 
skian is given by 

(2.8) W{6i{z),cj> l {z))=l. 

The eigenvalues of Hi with a Dirichlet boundary condition at x = 1 are given by 
the zeros of the entire function <pi{z 1 1) which are the squares of the positive zeros 
of the Bessel function J/ + i/ 2 of order I + 1/2: 

(2-9) Hl, n = (j'i + 1/2,™) 2 - 

Similarly, the eigenvalues of Hi with the boundary condition /'(l) + /3/(l) = 
are given by the zeros of the entire function </>J(z, 1) + (3(f>i(z, 1) and the positive 
eigenvalues are the squares of the positive zeros of zJi +3 / 2 ( z ) ~ (ft + 1 + l)Ji+i/2( z ) 
(cf. [HI §3.2]): 

(2-10) = (if + i/2,J 2 - 

Here the eigenvalues are counted according to \f < fi^i < < . . . . In particular, 
observe that the first eigenvalue Af will be zero for j3 = — (I + 1) and negative for 
0<-(l + l). 

Note that asymptotically ([TJ (9.5.12)]) 

(2.11) y/f*^ = 3i + $, n = ( n +£) ^ + 0(n~ 1 ), 

(2-12) v /^ = i f +|in =j J+iin + 0(n- 1 )= (n+^y+Oin- 1 ) 

for fixed I. 

Now let us look at perturbations 

(2.13) H = Hi+q{x) 

assuming that the potential q satisfies the following conditions: 

Hypothesis 2.1. Let I 6 [— 5,00). Set 

< 2 '"> '» = (n W !' mvmi ! >_ '' 

[(l-log(a;))|g(a;)| ) i = -£, 
and suppose q is real-valued such that 
(2.15) xg(ar) € £ x (0, 1). 
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Lemma 2.2. Assume Hvpothesis \2.1\ Then there is a solution tf>(z,x) of Hf = zf 
which is entire with respect to z and satisfies the integral equation 



(2.16) <f>(z,x) = <j>i(z,x)+ Gi(z,x,y)q(y)(j>(z,y)dy, 

Jo 

where 

(2.17) Gi(z,x,y) = <f>i(z,x)6i(z,y) - cf)i(z,y)6i(z,x) 

is the Green function of the initial value problem. Moreover, this solution satisfies 
the estimate 

(2.18, - < c ( TT ^)' +1 ^"■>"X T$fo>- 

The derivative is given by 

■h'Jy. t\ a- ( 

dx 



d 

(2.19) 4?(z,x) = <f/ l (z,x)+ I —G l {z,x 1 y)q{y)(t){z,y)dy 



and satisfies the estimate 



,2,0, |„ ( „„ - *(,„>, < c ( TT ± m -)'^> l , £ _iM_ d , 

Proof. In a way similar to [6] (see also |14j ) this can be shown by iteration using 
Lemmas I A. 1 1 and IA.2I Namely, it is easy to show that cj> given by 



(2.21) = 



H,n, 9l,0 ■= 91, 



n=0 



(2.22) (f>i t „ +1 (z,x) := / Gi(z,x,y)q(y)(j)^ n (z,y)dy, neN , 

Jo 

satisfies (|2.16[) . The inequalities 

« ^''^(TT^r^OTT**)"- 

which are necessary to prove the convergence in (|2.2ip . follow by induction from 
(|2TT5| and Lemma IA~T1 By (ET2Tj) and (pT2"5)) . we get inequality (j2~T5)) . where g is 
given in (I2.14[) . Hear C is to be understood as a generic constant whose value is 
different in (|2.18p and in Lemma I A. II The properties of <f>' follow in a quite similar 
way from (|2~2T|) - (|2~2"g|) and Lemma [O Finally, using ([23), (j2TT^l) . ([2~T7|) . and 
(|2.19p one can see that Hcj) = z<fi. □ 

Using the well-known asymptotic formulas for Bessel function [TJ (9.2.1), (9.2.11)] 
we obtain 

(2.24) tj>{z, x) = z-^ (sm(^z~x - y) + 0{\z\- l ' 2 e^ m ^)^ , 

(2.25) 4>'{z,x)=z-i (^cos^x-^+Odzr^VI 1111 ^)!)^, 
as z — > oo. Next, note that 



(2.26) lim x'^^z.x) = lim x~ l ~ 1 9i{z, x) 



11111 tXj \ls \ iO ■ tXj I 11111 Uij Vis I \ As « tAj I , j , -. • 

rv ' ; x^o YlK ' ; r(Z + |)2'+ 1 
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Moreover, using d'Alembert's formula (cf. [8j Sect. XI.6]) a second linearly inde- 
pendent solution, satisfying W(6(z), 4>{z)) — 1, is given by 

dy 



(2.27) 6(z,x) = -<j>(z,x) 



<l>(z,y) 



2 ' 



where c = c(z) has to be chosen such that <j)(z,x) does not vanish in (0, c). In 
particular it is straightforward to show 

Corollary 2.3. Assume Hypothesis \2.1\ The differential equation H f — zf has 
two linearly independent solutions, satisfying W(9(z),(j>(z)) = 1, of the form 

(, 28 , «,,„) _^,,«), t~-i 

where 4>{z, x), 9(z, x) £ C(C, [0, 1]) are jointly continuous and 4>{z, 0) = 9(z, 0) _1 ^ 
0. 

Unfortunately, since c(z) — > as z — > oo, this simple approach to get a second 
solution looses control over 9(z,x) as a function of z. In particular, it is not clear 
that this second solution can be chosen to be entire as a function of z. We will 
construct a second solution with better control with respect to z in Lemma | 
below. 

In any case, the behavior of the solutions of Hf = zf near x = implies 



Theorem 2.4. Assume Hypothesis \2.1\ The differential equation H = Hi + q is 
limit circle at x = if I £ [— i, |) and /imrf point at x — /or i > i. 7n particular, 
H associated with the boundary conditions at x = (/or ^ 6 [— i |) J and £ = 1 

(2.29) lW((Z + l)/(z)-z/'(z)) =0, /(l) = or /'(l) + 0/(1) = 

is self-adjoint. Moreover, the spectrum of H is purely discrete and bounded from 
below. 

Proof. Since (f)(z,x) and 9(z,x) are both in £ 2 (0, 1) if and only if / e [—5, |) we 
deduce that H is limit circle at x if and only if I € [— f> 2)- Moreover, in this case 
we can choose the boundary condition (cf. [TS1 Sect 9.2]) 

\im WM0,x),f{x)) = — lW((l + l)/(aO -*/'(*)) = 0, 

k->o I (/ + ^)2 1+i =c->-o 

where we have used (|2.26|) and the fact that any solution in the maximal domain 
of the differential expression satisfies lim x _j.o f(x) = 0. 

Furthermore, since (f>(z, x) has only a finite number of zeros inside (0, 1), the dif- 
ferential expression is nonoscillatory for every z G R, we conclude that the spectrum 
is purely discrete and bounded from below (cf. [El Thm. 14.9]). □ 

Now we are able to prove our eigenvalue asymptotics using a refined version of 
the approach by Poschel and Trubowitz |llj . 

Theorem 2.5. Assume Hvvothesis \2. 1\ Then the Dirichlet eigenvalues correspond- 
ing to the boundary condition /(l) = satisfy 

(2.30) fi n = (ji+^n + £„) 
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where 

Similarly the eigenvalues corresponding to the boundary condition /'(l) -1-/3/(1) = 
satisfy 

( 2 - 31 ) Ag= (if + , n + e ~„) 2 

where s n is of the same order as e n . 

y'n/2 z 4 J ;+1/2 

Then our estimate f|2 . 18[) reads 

e l Im V^I r 1 vd(v) 

\<Kz)-Mz)\<C- j==——e{z) , e(*):=/ - l™> dy. 

(i + y/\z\y +i Jo i + M 1 /^ 

Next, using [U (9.2.1), (9.2.11)], we have 



Proof. We set <f>i(z) := 0/(2,1) = 2 ^ Jz+i/ 2 (\/i) and </>(z) = 0(2,1). 



C>.:-!2) ,„(:): Un^z - j ) + 0(\z\~ V^IM* 1 ")! 

(2.33) = \*-^ (cos(V^- ^) + 0(|z|-V2 e IM^ 2 )l 

where the dot denotes a derivative with respect to z. Hence, taking into account 
(j2~TT]) . (j232|) . and ([233]) , we get 

(2.34) <M(^n) 2 ) = 0' ^(Oi + |,J a ) = ^ + §,J~ ( ' +2) ((- 1 )" + 0(0) . 

Furthermore, (|2.32p ~ (|2.34|) together with the mean value theorem yield (for suffi- 
ciently large n) 

IMUl+hn ± £ ») 2 )l ^ \tit+$,n + en)- il+2) \(j l+hn ± £„) 2 - jf+ij 

. ien|2j I+ i in -e n | l e „ e C/i+*,n) 
> ttt 2 rrrr > tt: rrrr = 2(7- 



where e„ = 8Ce(j ;+ i „). Thus 

\m +h n±Cn) 2 )-MUl +h n±Cn) 2 )\ <C '^Z^L 

< |^((i /+ i„±£„) 2 )|. 

This shows that (f>(z) has different signs at (j l+ i_ n — e n ) 2 and (j ;+ i.„ + s n ) 2 , and 
thus there is at least one zero in between. 

Let us show that </> has no other zeros. Since | sinz| > ie' 1 ™ 2 ' if \z — nn\ > j for 
all n £ Z (see for instance [TTJ Lemma 2.1]), there exists a N £ N such that 



e' 1 ™^! lir , 7r 

Z 2 ^ 



with some positive constant Ci > independent of 2 (this estimate can also be 
deduced from [9J Lemma 22.1]). Since e(z) — > as |z| — > oo, there exists K > 
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such that e(z) < (C;C) _1 if \z\ > K. Thus, on contours \z\ = (K + N + 1 -4^)tt and 
\y/z~ rrn - l -f\ = j, n> K + N +1, we obtain 

\4>{z) - 4>i{z)\ < \<fn{z)\. 

By Rouche's theorem, 4>{z) has as many roots as 4>i(z) in each of the bounded 
regions and the remaining unbounded region. Since all roots are simple, we are 
done. 

An analogous argument can be given for the general eigenvalues based on zeros 

of 



<pf{z) = ti(z, i) + i) = \j\z-^ ((0 + 1+ 1) J l+h - <JZJi +i (Vzj) ■ 

□ 

Similar results for q £ L 2 were given by Guillot and Ralston [6 , Carlson [4], 
Serier [14 . Note that the case xq(x) £ W for I £ N is covered in [2 E2 H3] for Z = 
and in [3] for Z £ No, where the direct and inverse spectral problems have been 
effectively studied for operators with distributional potentials q £ W _1 ' p (0, 1). 

Finally we come to uniqueness results for the inverse problem. Considering the 
solution 

(2.35) ^{z,x) = e l {z,l)<l> l {z,x)-<ln{zA)9i{z,x)=Gi{z ) x,l) 
satisfying the initial conditions (tj}f (z, 1), (ip^°(z, 1))') = (0, 1) and 

itf(z,x) = -(m(z, 1) + e[{z, l))<fn(z, x) + (pMz, 1) + 4>[{z, l))6i{z,x) 

BC 

(2.36) =-pG l (z,x,l) + —±(z,x,l) 

dy 

satisfying the initial conditions (ipi(z, 1), (ipi(z, 1))') = (1, — /?), we obtain the ana- 
log of Lemma 12.21 

Lemma 2.6. Assume Hypothesis ] 2. 1\ Then there is a solution ip@(z, x) of Hu = zu 
which is entire with respect to z and satisfies the integral equation 

(2.37) ^(z,x)=^(z,x)- ( G l (z,x,y)q(y)^(z,y)dy. 

Moreover, for I > —1/2 this solution satisfies the estimate 
(2.38) 

/l + \z\ 1 / 2 x\ l e l Im 1/2 )l( 1 - iB ) r 1 ua(u) 

( , 39 , ^ I) . rf( ,, I) | SC (^)' J M^<-. J r_^ 5 .*. 

In the case I = — ^ an additional factor (1 — log(a;)) appears in the right-hand sides 
of (I2~38l) and (£Q9l) . 
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Proof. Suppose /3 oo first. In a way quite similar to the proof of Lemma 12.21 one 
can show that the solution ipP admits representation 



(2.40) ^ = X>L <o-=W> 



n=0 



(2.41) i>l n+1 {z,x) :=- / Gi(z,x,y)q(y)tf n (z,y)dy (n G N ), 

where the functions satisfy inequalities 

for Z > — K and inequalities 

(2.43, | <i( «>|<^ (1 _ log(l)) (£±^) ,/2 

Vix l + M^y 

for / = — i. Indeed, the inequality (|2.42j) for n = and for some C > follows 
from (|2.36j) . (|A.22p . and the fact that according to (|2.17p we have Gi(z,x,y) — 
—Gi(z,y,x), Next, we use (|2.41j) and (|A.2I) to prove (|2.42l) for all n by induction. 
By (|^3D|) - (|02]) we get ([237)1 and (|2~M1) . Finally, according to (|2~37|1 we have 
= z^ and (^(z, 1), (^(z, 1))') = {^{z, 1), (Vf (*, 1))') = (1, -/?). 

The inequality (gj5]) for n = follows from (|2~M1) and (fA~23l . and the remaining 
part of the proof for I = —1/2 is analogous to the case I > —1/2. 

The case j3 = oo can be treated in a similar way, it suffices to note that 

(2 .44, <_^(i±^!£)'e,M,-,,<^ 
for / > — i and 

(2.45) m^ Z , X )\ < (l-lo g (,) )TT ^ (^|^) 1/2 e |Im(zl/2)l(1 ^ 
for I = i. □ 
Remark 2.7. 7n [18] a second solution ip(z,x) is constructed by considering 

(2.46) ^(z,x) = 9i(z,x) - Gt(z,x,y)q(y)<tp(z,y)dy. 

J X 

It is claimed to be linearly independent for all z and a reference is made to [6], 
where the corresponding claim for I — 1 was made. However, this is wrong as the 
following counter example shows (a similar counter example was communicated to 
us by Ralston; see also ^7\): 

The function ip satisfies (ip{z, l),?//(;z, 1)) = (0i(z, l),9'i(z, 1)). Let yi+1/2,11 be a 
zero of z) and chose q(x) such that yi+i/2, n * s an Dirichlet eigenvalue of H. 
Then W(tp(yi +1 / 2 , n ), ^(y/+l/2,n)) = 0. 

Now we come to our uniqueness results: 
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Theorem 2.8. Assume Hypothesis \2.1\ The following set of spectral data deter- 
mine q and I uniquely: 

(i) Two sets of eigenvalues A" and A^ for a ^ (3 (including the case A^° = [i n ). 

(ii) The Dirichlet eigenvalues {/in}^i together with the norming constants 

1 f 1 At.. „\2r_ HVnA) 



7n = T77 T\T / 4>{Hn,x) dx 



0'(>ri,l) 2 Jo ' 4>'{P-nA)' 

or the set of eigenvalues {A^}^L together with the norming constants 

where c^(z) := ^'(z, 1) + /3</.(z, 1). 
(iii) The eigenvalues A^ together with </>(A^, 1) if (3 ^ oo or <fr'(\P, 1) if ^ 0. 

Here the dot and prime denote derivatives with respect to z and x, respectively. 

Proof, (i) To see the first claim recall the Weyl m-function of H at the regular 
endpoint x = 1, which is given by 

N _ <Kz,l)-/V(z,l) 
mM2) -0'(z,l) + /3^,l) 

and consider 

, , N a/3 + 1 l + (3 2 (j)(z,l) + a^(z,l) 

771 Z — 771/5 Z) = • ; : ; r. 

( ' p{ ' (3 -a a-0 (P'(z,l)+I3<p{z,l) 
Then the well-known asymptotics [15] 



TO/3 (z) 



- v / ^+0(l), /3 = oo, 

+ o(i), /3eR, 



show that rh(z) is uniquely determined by its zeros and poles via its Hadamard 
product. To complete the proof of (i) it suffices to note that the Weyl function 
mp(z) uniquely determines the potential q and the constants I and {3 (see |10jV 

(ii) Similarly, the Herglotz function mp(z) is uniquely determined by its poles 
A^ and residues —7^. 

(iii) Let (3 ^ 00 and 0(A£, 1) be given. Clearly I can be read off by the asymptotics 
of \@. So let us consider H and H associated with q and q, but with the same I. 
Now, following [111 Chapter 3] consider the meromorphic function 

(<j)(z, x) - 4>(z, x)) (tpPjz, x) - TpPjz, x)) 
The only poles are at the eigenvalues A^ and the residues are given by 

(^\i,x)-k\ p n ,x)y 



, 2 

>o, 
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where we have used (j)(X^,x) — <p(X^, 1)^ /3 (A^, x) together with our assumption 
<f>(\%, 1) = 4>(\P, 1). Using our estimates (j2T8| . pT20|) . and (TT39"]) . we have 

(*,*) - *(*,*)! < C [ ) m e IM^)|. f ^+P ^ ; 



|<c( 




a; 




vl + 


M 


l / 2 x) 


|<c( 


'1 + 


M 




V X + 




iVaaj J 



l + lzl 1 ^ 

- ^(z,*)| < C ( l±l£l^) l e IM^)l(^) Z 1 

Further, together with the estimate from the proof of Theorem 12.51 we get 

e IM*) 1/2 l 



\^{z)\ > \tf{z)\ - \<f(z) - tf{z)\ > \\tf{z)\ > Ci- 



2 2 



for \z\ = r n = (n + (1 + l)/2) 7r 2 with n sufficiently large. Observing that for every 
x e (0,1) 

1 V fl\ + l y)) dy < f + ^ = °d^ 1/2 ) 

l + \z\ l -/ 2 y l + \zy-i z x J x 

we see that the function (|2.47l) is o(r~ 1 ) along these circles and thus Lemma 3.2 
from shows that the residues are all zero. Hence <fi(\f,x) = 4>(X^,x) implying 
q = q. 

If (3 and the set </>'(A£, 1) is given, then the proof remains essentially the same. 
The only difference is that we now use </>(A^, x) — —(3^ 1 (f>'(X^, l)i/^(A^, x) if (3 ^ oo. 
In the case (3 = oo we need to use <p^(z) = <f/(z, 1), <p(X^,x) — (f>'(X^, (X^,x), 
and the estimate (f2T38|) instead of (f!Q9]) . □ 

Appendix A. Some estimates for the spherical Schrodinger equation 

In this appendix we want to provide some estimates for the solutions of the 
spherical Schrodinger equation which are crucial for the main body of our paper. 
These results are due to Guillot and Ralston 6 in the case 1 = 1. The analog 
estimates for arbitrary I > — | have been stated in [18] without proof. However, 
the estimate (39) from [18] is clearly wrong in the case I = — |. The case of integer 
I is given in [14]. Since the proof in [14] (and [6]) uses the explicit representation of 
spherical Bessel functions in terms of trigonometric functions we have decided to 
provide the details for the general case in this appendix. We assume that x <G (0, 1). 

Lemma A.l. For I > — i the following estimates hold: 

(A.i) l^,*)l<c( TT ^—f^K 

(A.2) |«K^,ril<c( T ^)' +I (l±J^)'el-<.-»K.-.», y < x . 
For the case I = — i formula (jA.l[) remains valid and one has to replace (|A.2I) by 

(A.3) |G_ 1/2 (z, x, y)\ <C [ {l + + ) 

X e IM* 1/a )l(*-»)( 1 _ log(y)), y < x. 
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Proof. First of all recall fljO]), flU}, and (|27T7|) . and note that for all I > -\ 

(A.4) <f>i(z,x) = A^z-^J l+ ^{yfzx), 

G t (z,x,y) = --^y/xy{Ji+^{Vzx)Y l+ i(s/zy) - J l+ i{J~zy)Y l+ i(Jzx)) 
(A.5) = -^(^i {V~zx)H^ h (A~zy) _ Hj W (Viy)^ (Via:)) . 

Here ff, ( + 1/2 (z) = Jj+i/aC^) +1^+1/2(2) and #4+1/2 ( z ) = J ;+i/2( 2 ) - 1^+1/2(2) arc 
the Hankel functions of the first and second kind, respectively. Moreover, recall 
Yi +1/2 (z) = sin((2 + l/2)^)- 1 (cos((Z + 1/2)tt) J l+1/2 (z) - J-i-i/ 2 (z)) (ffl (9.1-2)]). 

The claimed estimates (|A.1I) and (IA.2I) can be shown by combining the following 
asymptotic expansions which yield estimates for \z\ > 1 and \z\ < 1, respectively. 
In fact, we have [H (9.2.7)-(9.2.10)] 

(A.6) H^Az) = J— e i ( 2 -('+ 1 )' r / 2 )(l + 0(l/z)), -tt < arg(z) < 2tt, 

'+2 V 7TZ V ' 

(A.7) ffS(z) = J— e-K^-C+i)^) h + o(l/z)), -2n < arg(z) < tt, 

'+2 V 

for z — > 00 with the error uniform strictly inside the indicated sectors ( |171 (7.2.1)- 
(7.2.2)]). Using the power series for the Bessel and Neumann functions Q~l (9.1.10)- 
(9.1.11)] we get 

(A.8) J l+k {z) = (|) /+ " (1 + O(z)), 

(A . 9) r-^^)--* a+0 (^% ;+ i>o, 

[|log(2) + 0(l), I = -3« 

for z ->• 0. Thus, (|A.1|) immediately follows from (|A.6|) . (|A.7|) . and (IA.8I) . The 
estimate (IA.2I) can be also deduced from (|A.6[) - (IA.9I) . For instance, using (IA.8I) . 
(|A.9|) . and the fact that the function x h+ x(l + x)^ 1 is increasing, we get for 

m < lei < 1 



< c 




lei 



Similarly one handles the cases \rf\ < 1 < |£| and 1 < \rj\ < |£| to obtain the desired 
result (for the last one use the expression in terms of Hankel functions) . 

The estimate (|A.3j) requires some further considerations, and we also split this 
case into three subcases: a) sufficiently large values of \z\ 1 / 2 y] b) the values of 
\z\ x / 2 x are sufficiently large and the values of | 1 1//2 2/ are bounded; c) \z\^l 2 x is 
bounded. In all the subcases it is assumed that y < x. 
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To prove (jA.3|) for the subcase of sufficiently large values of Izl 1 / 2 ?/ one can use 
again the equality (|A.5|) : 

(A.10) G_ 1/2 (z,x,y) = ~^xy(H^(^x)H^(^y) ~ H^i^H^^x)). 

Taking into account (| A.6|) and (|A.7|) we easily derive (|A.3j) in the subcase a). 
Assume now that the values of \z\ 1 ^ 2 y are bounded. Note that (see (|2.5p . (|2.17j) ) 

G-i/ 2 (z,x,y) = 4>_i(z,x)9_i(z,y) - <p_i(z,y)6_i(z,x), 

where 

9_i(z,x) = (y q {s/zx) - ^log(z)J (V^a;)^ . 

Using the integral representations [TJ (9.1.18) and (9.1.19)], we have 

(A.ll) 6_i{z,x) = ir- 3/2 V2x~ [ cos (y^x cos d) {-y + \og(x) +\og(2 sin 2 9)} d6. 
2 Jo 

It is immediate from (|A.11|) that 

(A.12) \6_i(z,x)\ <dV^{l~log(x))e llmizl/2)lx . 

Further, since the values of |z| 1//2 i/ are bounded, by (|A. 1 1) and (|A.12[) we derive 

/ \ 1//2 

(A.13) 1 (z, x)6_ k (z, y)\<C 2 ( TT ^) (1 - log( Z/ ))el Im ^ 1/2 )l- 

MjVTWW^T* < "" ,l< - > ('-^))- 

For large values of |2r| 1,/2 a:, Hankel's asymptotic expansions [H (9.2.5) and (9.2.6)] 
yield: 

(A.14) \6_i(z,x)\ ^^Izl-^^l+logdz^el 1 " 1 ^ 172 '!". 

Since \z\ x l 2 x is sufficiently large, the values \z\ and log(|z|) are also large. Moreover, 
\z\ x / 2 y is bounded and hence log(|z|) + 2 log(y) is bounded from above. Thus, ljA.14j) 
can be rewritten in the form 

(A.15) |«L*(*,*)| < C 5 f — _ j e l Im ^ 1/2 )l-(l -log(y)). 

It follows from (|A.1|) and (|A.15|) that 

\ 1/2 

(A.16) \8_i(z,x)0_i(z,y)\<C e < 



(l + |z| 1 /2 a; )(i + | z |i/2y) 

x el Im ^ 1/2 ^ x -^(l-log(2/)). 

Inequalities (|A.13I) and (|A. 16|) yield (|A.3[) for the subcase b). 

Finally, if Izl 1 / 2 ^ is bounded, according to (|A.1[) and (| A. 12|) we have 

(A.17) \e_x_{z,x)^{z,y)\<C 7y fxy{l-\og{x))<C 7y /xy{l-\og{y)). 

Inequalities (|A.13I) and (|A.17|) imply (IA.3I) for the subcase c). □ 
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To consider J^0z(z,x) and -g-Gi(z,x) one should simply use the asymptotic 
relations [TJ (9.2.11)-(9.2.16)] on the derivatives of the special functions, in addition 
to the relations [TJ (9.2.5)-(9.2.10)] which were used in the proof of the previous 
lemma. In particular, for sufficiently large values of \z\ we have 



(A.18) 

and, as z — > oo, we get 
8 



9 J 



< Cs\ — e^WI, 



(A.19) i-H}]\(z) = iJ— e»(*-(l+i)*/2) (1 + 0(1/2)), -7T < arg(z) < 2tt, 
(A.20) ^-H^\(z) = -id — e- i(2 - ((+1) " /2) (l + 0(l/z)), 



7TZ 

Lemma A. 2. For I > — \ the following estimates hold, 
d 



~2ir < arg(z) < ir. 



(A.21) 
(A.22) 



dx 



4>i(z,x) 



-^G t (z,x,y) 



< C 



< c 



1 + \z\l/2 X 



JIm( 2 1 /2 ) | x 





x \ 




-N| 1/2 2/\ 


(n 


- \»\V*x) 







y <x. 



For I = —1/2 formula (IA.21[) remains true and one has to replace (|A.22I) &?/ 
(A.23) 



-^-G_ 1/2 (z,x,y) 
ox ' 



< 



y <x. 



Proof. The proof is similar to the proof of Lemma I A . 1 1 and we shall prove here only 
formula (|A.23|i . For the subcase of sufficiently large values of \z\ x / 2 y we use the 
equality (|A.10I) . Therefore, in view of (IA.6I) . (|A.7j) . (|A.19|) . and (|A.20|) we derive 

(A.24) |G_ 1/2 (z,x, 2 /)|<C7 9 el Im ( zl/2 )l^), y < x, 

and (j A.23|) is immediate. 

When the values of Iz^^y are bounded, formulas (|A.12[) and (| A.21I) imply the 
inequality 

(A.25) 

1/2 



d 

—<f>_i(z,x)9_i(z 7 y) 



< 



C W ( V ^Z XV ) elM^)l(^) ( i_ log(y)) , 
\X -\- \z\ xy J 



which we shall use instead of (|A.13[) . For the subcase of sufficiently large values of 
\z\ 1/2 x and bounded values of |z| 1/2 ?/, it follows from (j2~5j) . (|A.15|) . and (|A.18I) that 

|l/2_\ V2 



(A.26) 



_9_ 



; (z,x) 



<di 



1 



3 l I ^ I/3 )l»(l-log( 2/ )). 



Taking into account (|A.1|) . ()A.26|) . and boundcdness of \z\ x ^ 2 y we have 
(A.27) 



d 

—6_i(z,x)<j)_i(z,y) 
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Inequalities (|A.25j) and (|A.27|) yield (|A.23j) for the subcase of sufficiently large 
values of \z\^l 2 x and bounded values of \z\ 1 ^ 2 y. 

Finally, if the values of \z\^-l 2 x are bounded, formula (|A.11|) implies 



(A.28) 



dx 



Az,x) 



^Cux'^il-logix)), 



and (|A.27j) follows from (|A.1|) and (|A.28j) . By (|A.25j) and (|A.27j) the inequality 
(|A~23)) holds again. □ 
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